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Introduction 

The  Laplace  inversions  treated  herein  arise  in  connection  with  some 
conjugated  heat  transfer  problems,  details  of  which  have  been  reported  in 
a  thesis  by  one  of  the  authors  [4]  and  will  also  be  reported  elsewhere, 
as  well.  In  connection  with  another  physical  problem,  R.  A.  Kantola  [l] 
arrived  at  similar  results.  However,  the  details  here  are  sufficiently 
involved  and  sufficiently  different  than  those  dealt  with  by  Kantola  to 
suggest  preserving  our  analysis  in  a  monograph  such  as  this. 

In  the  body  of  this  monograph,  a  theorem  is  stated  and  proved  which 
affords  the  Laplace  inversion  of  a  rather  general  function.  Appendix  A 
indicates  a  suggested  numerical  procedure  for  evaluation  of  the  definite 
integrals  appearing  in  the  result.  Appendix  B  treats  certain  details  of 
a  special  but  useful  case.  Appendix  C  merely  lists  some  additional  Laplace 
inversions  which  have  come  to  our  attention  in  connection  with  related  work 


but  which  do  not  result  from  the  analysis  given  in  the  body  hereof. 

Statement  of  Theorem 
The  principal  content  of  this  monograph  is  contained  in  the 
THEOREM:  Let 

f(s)  =  \  exPL--/-ifsT  ]  (1) 

where  m(s)  is  a  meromorphic  function  whose  zeros  s  =  -a.  and  whose  poles 


s  =  -b.  satisfy 


o  <  a-,  <  b-,  <  a2  <  b?  <  . . .  <  aN  <  b».  ^  oo        (2) 


Then,  for  t  >  0,  the  inverse  Laplace  transform  F(t)  =«0  (f(s)} 


is  given  by  N  b. 

« 1  1       K 


F(t)   =  exp  [-y(o)]   "  7    )  /         7exP   (-xt)   sin  [M(x)]  dx         (3) 

where  k 

u(x)  =  /|m(-x)l  (4) 

The  sum,  having  upper  limit  N,   is  to  be  interpreted  as  follows. 
If  there  is  a  finite  number  N  of  pairs  a,    <  b, ,  no  special   interpretation 
is  necessary.     If  there  is  a  finite  number  (N-l)  of  pairs  a.    <  b.    and  an 
additional   zero  s  =  -a,,,  a..  >  bN_, ,  an  Nth  pole  s  =  -bN  =  -  °°  is  to  be 
added  so  as  to  make  N  pairs.     If  there  is  an  infinite  number  of  pairs 
a.    <  b,  take  N  =  °°. 

Proof  of  Theorem 
In  the  proof  we  consider  a  sequence  {R  }  of  positive  numbers, 
R    t  &\,>  b.    for  any  n  and  any  k,  R    ->  »  as  n  ■>  °°,  and  we  consider  a 
corresponding  sequence  of  closed  contours  C     in  the  complex  s-plane  each 
of  which  encloses  the  pole  of  f(s)  at  s  =  0  but  excludes  all   of  the  branch 


points  at  s  =  -a.  and  s  =  -b..  The  contour  C  consists  of  four  sub- 
contours  C  ■.,  C  p»  C  o»  anc'  ^n4  as  follows.  C  ,  is  a  vertical  line  from 
s  =  -i  R  to  s  =  i  R  ,  indented  as  shown  by  a  small  semicircle  passing 


Figure 
Contour 


to  the  right  of  s  =  0.  Cn2  is  the  arc  s  =  Rne10,  f  =  e  £  *  and  Cn4  is 
the  similar  arc  s  =  Rneie;  -n  i   e  i   -  J.  Finally,  C^  is  the  contour 
which  passes  to  the  right  from  s  =  -Rn,  along  the  negative  real  s-axis 
to  s  =  -a,  with  small  semicircular  indentations  into  the  upper  half 


s-plane  about  each  zero   (-ak)  and  about  each  pole  (-b, )  after  which  it 

then  passes  to  the  left  from  s  =  -an    to  s  =  -R     with  small   semicircular 
r  I  n 

indentations  into  the  lower  half  s-plane,  as  indicated  in  Figure  1. 


Clearly,  by  Cauchy's  theorem 

r 


Jn  =  7-5-     i  eSt  f  (s)  <*s  •  <WP  [-/TSTOT  ]  =  exp  [-u(o)]       (5) 


1 
2iri 

for  n  =  1 ,   2,  3,    . . .  Also,   by  the  complex  inversion  theorem  for  Laplace 


transforms  ^ 

Jnl   =  ST    \       eSt  f(s)  ds  *F(t)  (6) 

as  n  -*  °°. 

To  proceed  further,  we  wish  to  establish  an  inequality  for  arg[m(s)]. 
To  do  this,  we  consider  complex  vectors  emanating  from  the  points  s  =  -a. 
and  s  =  -b.  ,  and  we  write 

1    +  s/ak  =  uk  exp(i'4>k);  uk  >  o;   "*  =  \  =  *  (7) 

1   +  s/bk  =  vk  exp(i>k);  vk  I  q;   -tt  =  ^k  I  tt  (8) 

(s)  =  m(o)  n(n)   (1   +  s/aji  n(n)   (1   +  s/b.)  (9) 


mn 


where 


and 


u 


=  ^  (s)  exp  (ien)  (10) 


2  (s)  =y(0)  Ji{n)  uk  *n<n>  vk  (11) 


The  sums  and  products  written  with  a  superscript  (n)  indicate  that 

only  those  terms  are  to  be  considered  for  which  a„  <  R    or  b     <  R  ,  as  may 

J  n         n  n         n 

be  appropriate. 


Now  note  for  any  s  in  the  (closed)  upper  half  plane  we  have 

4>k  I     i^,       k  =  1,   2,    ...  (13) 

A  geometrical  argument  is  given  in  Figure  2.  Now  either  there  are  N 

pairs  a.  ,  b,  each  less  than  R  ,  so  that 

N 

en=   X>k-*k)>=0  (14) 

or  there  is  an  extra,  unpaired  zero,  so  that 

N-l 
en=   y>k-\>+*N  =  0  05) 


s  =  s 


Figure  2. 
Demonstration  that 

♦k>*k 


In  the  first  case,  we  can  write 

N-l 

Y,  (*  k  "  *k+l }  '  ^N  "  *1   "  *N  <  *1   =  ff 


'n       *! 


k=l 


whereas  in  the  second  case  we  can  write 

N 

en  =  *1  _   E  (*k  "  W  =  *1  =  * 


Thus,  in  either  case, 


k=l 


0  =  6n  =  * 


(16) 


(17) 


(18) 


Similarly,  by  considering  an  s  in  the  lower  half  plane,  we  find 


-.  S  6n  i    0 


(19) 


so  that  in  all   cases 

-TT    =     3N    i    TT  (20) 

Thus,   if  we  write  r 

\l"h      \        estf(s)ds  (21) 

Cn2 
we  note  that 


f(s)    |   =  J-  |   exp  (-  /  u  2  Tb       ) 
n 


^-exp  [-un  cos(Bn/2)]  l  ^~ 


(22) 


uniformly  on  C  «»  °y  virtue  of  the  inequality  (20).     Thus,  by  Jordan's 

lemma  J  „  +0  as  n-**>.  Similarly,  the  integral   on  C  ,  vanishes  as  n  ->  ». 

Next  we  consider  the  integrals  around  semicircles  surrounding 

poles  and  zeros  on  contour  C  ~.     The  result  is  obvious  near  a  zero: 

s  =  -a,;  es     -  e"ak  ,       f(s)       -  —    since  m(s)  -  0,  and  the  length 
k  ak 

of  the  arc  goes  to  zero  as  its  radius  does.     Thus,   the  integral   on 
such  a  semicircle  about  a  zero,  s  =  -a.  ,  vanishes  with  the  radius. 
However,  for  the  integral   on  a  semicircle  about  a  pole,  s  =  -b,  ,   it 
becomes  necessary  to  use  a  considerably  more  subtle  argument.     For 
this  purpose  we  must  invoke  (20)  again.     We  have 


I   exp[-m(s)]    |    =  exp[-Mn  cos(en/2)]  <  1  (23) 

even  though  y     ->  °°  as  the  radius  of  the  indentation  goes  to  zero. 
Thus,   the  integral   on  such  a  semicircle  about  a  pole,  s  =  -b. , 
vanishes  with  the  radius. 


We  remark  that  although  we  have  reason  to  believe,  on  physical 
grounds,  the  correctness  of  an  evaluation  by  Kantola  [l],  we  feel 
that  his  argument  intended  to  show  vanishing  of  such  an  integral  on  a 
semicircle  surrounding  a  branch  point  is  deficient.  Our  argument, 
based  on  (20),  is  not  applicable  to  his  case. 

If  now  we  write  (5)  in  the  form 

Jnl  "  Jn  "  Jn2  "  Jn3  "Jn4  <24> 

and  let  n  ->  °°,  and  also  let  the  radii  of  all  semicircular  indenta- 
tions on  C  3  go  to  zero,  we  find 

F(t)   =  exp[-y(0)]  -  J3  (25) 

where 

J3  =     2^~ 


JU-      /       jexp[st-v¥TsT~  ]  ds 

J  C 
L3 


and  C3  now  is  a  contour  extending  from  s  =  -~  along  the  upper  edge  of 
a  branch  cut,   to  the  right,  passing  around  s  =  -a- ,  and  then  extending 
along  the  lower  edge  of  the  branch  cut  to  s  =  -°°.     The  branch  cut  is 
the  negative  real    s-axis  from  -°°  to  -a,.     There  are  no  longer  any  semi- 
circular indentations  in  C~. 

To  reduce  (26)   to  a  more  convenient  form,  we  first  introduce 
the  function  v(x),  for  real   x,  which  denotes  the  excess  of  the  number 
of  zeros   (s  =  -a.  )  over  the  number  of  poles   (s  =  -b,  )  to  the  right  of 
the  point  s  =  x.     It  is  clear  that  v  =  0  in  the  interval   -cIl+i   <  s  <  -b. 
while  v  =  1    in  the  interval   -b.    <  s  <  -a,  .     Writing  6  now,  rather  than 

6   ,  since  we  have  taken  n  -»-  °°  and  all   zeros  and  poles  are  included  in 
n 

our  accounting,  we  note  that  on  the  upper  edge  of  the  cut  B  =  W2 


while  on  the  lower  edge  3  =-vtt/2.  We  also  now  write  y  rather  than  y  . 

n 

Thus, 

1       r\ 

J3  =  2^T  /  J  exp[st  -  y  exp(iv7r/2)]ds 

—00 

+  27T  s  exp[st  "  y  exP(-iW2)]ds  (27) 

-a] 

=     2~f  /  ?  eS      I    exp^"y  exp(ivfr/2)]  -exp[-y  exp(-ivTr/2)]  |    ds 

—00  V  -' 

If  v  =  0,   the  quantity  in  braces  vanishes,  while  if  v  =  1,   it  becomes 

exp(-iy)   -  exp(iy)  =  -2i   sin  y  (28) 


Thus 


-a 


k 


J3  =  -  7     2_        /  i"eSt  Sin  y  dS  (29) 

k-1       J-K 

If  we  now  change  the  variable  of  integration  from  s  to  x  =  -s, 

we  obtain 

N  b 

/        I 
[e-xt  sin[y(x)]dx  (30) 

and,  thus,  finally 

v-  rbk 

F(t)  =  exp[-u(o)]  -  \  )  /  I  e"xt  sin[p(x);|dx 

The  use  of  the  upper  limit  of  summation  is  consistent  with  the  remarks 
following  the  statement  of  the  theorem,  the  proof  of  which  is  now  complete. 
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APPENDIX  A 


Numerical  Evaluation  of  the  Integrals 


The  integrals  appearing  in  the  inversion  theorem  (3)  are  of  the 


form 


where 


I  =        h(x)sin[g(x)]  dx 
J  a 


(A-l) 

(A-2) 
(A-3) 
(A-4) 
(A-5) 


h(x)   is  continuous  in  [a,b] 

g(x)   is  monotonic  continuous  in  [a,b) 

9(a)  =  0 

g(x)  ->  °°  as  x  -+  b" 
It  is  contemplated  that  b  may  be  finite  or  infinite. 

The  integrand  is  illustrated  in 
the  figure  at  the  right.  The  integral, 
I,  is  the  algebraic  sum  of  the  shaded 
areas,  with  signs  taken  as  indicated. 
For  values  of  x  sufficiently  close  to 
b,  these  contributions  will  alternate 
in  sign  so  that  I  may  be  regarded  as 
the  sum  of  an  alternating  series  of 
terms  having  decreasing  magnitude. 
Thus,  the  error  of  truncation  is  less  than  the  magnitude  of  the  first  term 
to  be  neglected. 

In  the  numerical  treatment  one  must  evaluate,  in  succession,  the 
numbers  C  such  that 


Figure   A-l 


g  (Cn)   =  nir 
The  values  C     evidently  satisfy 


a  =  C~  <  C,   <  Cp  <   .    .    .   <b 


(A-6) 


(A-7) 
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This  determination  is  uniquely  possible  because  the  monotonicity 
of  g(x)  assures  the  existence  of  a  unique  inverse  function.     Many  tech- 
niques are  available  for  this  evaluation.     Some  remarks  suggesting  a 
procedure  suited  to  the  problems  which  form  the  subject  of  this  mono- 
graph will  appear  near  the  end  of  the  present  appendix. 

Thus,  we  have  N 


where 


I   =  Lim    V      qk  (A-8) 

C, 


qk  = J    h(x)  sin[g(x)]  dx  (A-9) 


'k-1 
Numerical   quadrature  is  indicated.     A  formula  of  Gaussian  type 

might  be  appropriate  except  that  the  ordinary  Gaussian  points  are  all 
interior  to  the  interval  of  integration  and  thus  cannot  reflect  our 
knowledge  that  sin[g(C  )]  =  0  for  n  =  0,1,2,...     A  slightly  more  ef- 
ficient procedure,  accordingly,  employs  Gauss-Radau  integration  in 
which  the  end  points  of  the  interval   are  specifically  included  among 
the  integration  points.     The  integrand  is  zero  at  these  points  so  that 
the  corresponding  summands  may  be  omitted.     For  details  of  Gauss-Radau 

formulas,   see  Z.    Kopal    [2].     The  appropriate  formula   is 

1  N-l 

F(x)dx  =  H^C-l)  +   £  H.   F(a.)  +  Hn  F(l)  (A-10) 

-1  J'=z 

where  a.  =  -a  _.+1  and  H.  =  H  -+,  are  tabulated  in  Appendix  4.5  of 

Kopal  [2]  for  N  =  2, 3,..., 11.  (Note  that  these  a.'s  should  not  be  con- 

fused  with  a.  's  which  locate  the  zeros  of  the  function  m(s)  employed  in 

the  body  of  this  monograph). 
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For  example,  for  n  =  5,  we  have 


GAUSS- 

■RADAU 

CONSTANTS 

FOR  N  =  5 

a5  =  1    =  -3] 

H5  =  H]    =  1/10 
H4  =  H2  =  49/90 

a4  =  /  3/7     =  -a2 

a3  =  0 

H3  =  64/90 

We  have  calculated  the  coefficients  a.  and  H.  for  the  case  N  =  12 

J     J 

and  list  the  values  here  so  as  to  supplement  the  information  given  in  [2] 


GAUSS-RADAU  CONSTANTS  FOR  N  =  12 

i]2  =  1.00000  00000  00000  0  H]2  =  0.01515  15151  51515  2 

in    =  0.94489  92722  22882  0  H,.,  =  0.09168  45174  13204  7 

i1Q  =  0.81927  93216  44007  4  H1Q  =  0.15797  47055  64367  1 

=  0.63287  61530  31861  8  HQ  =  0.21250  84177  61024  6 

9  9 

8     =  0.39953  09409  65348  9  HQ  =  0.25127  56031  99201  1 

iy     =  0.13655  29328  54927  6  H?  =  0.27140  52409  10696  2 


These  formulas  lead  to  exact  evaluation  for  a  polynomial   F(x)  of 
degree  2N-3.     In  our  case,  with  appropriate  mapping  of  [C   ,  C     ,]  onto 
[-1 ,  1],  we  have  F(-l)  =  F(l)  =  0  so  that  using  N  =  9,  say,  requires 
seven  nontrivial   evaluations,  leading  to  exact  integration  of  a  15-degree 
polynomial,  whereas  with  ordinary  Gaussian  integration,  we  would  have  had 
to  use  N  =  8,  i.e.,  eight  nontrivial   evaluations  for  exact  integration  of 
a  15-degree  polynomial. 
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The  actual  evaluation  of  q,  thus  becomes 

N-l 
qk  =  \   (Ck  -  CM)  )   Hj  •  h(Xj)  sin  [q(Xj)]    (A-ll ) 

j«2 


where 

xj4[<ck  +  ck-i)  +  dj(ck-ck-i)j  <A-12> 

The  value  of  N  to  be  used  should  be  consistent  with  the  accuracy  re- 
quired and  the  costs  of  computation 

To  test  this  procedure,  we  considered  the  special  case  where 

h(x)  =  (1/tt-x)"2  exp[-l/(l/Tr-x)]  (A-13) 

q(x)   =  ttx/O/tt-x)]  (A-14) 

a  =  0;  b  =  1/tt  (A-15) 


We  easily  find 


C     =  n/^(n+l)  (A-16) 


and  using  N  =  6,  we  obtain  the  following  evaluation 

*  ck  qk  Eqk 

0  0- 

1  1/2tt  .0225407  .0225407 

2  2/3tt         -.0009741  .0215666 

3  3/4tt  .0000421  .0216087 

4  4/5rr         -.0000018  .0216069 

5  5/6tt  .0000001  .0216070 

In  this  particular  case,  the  integral  can  be  evaluated  exactly.     We 
have  I/* 

I  =  tt2     /  (1-^x)"2  exp[-Tr/(l-7Tx)]  sin  [>2x/(l-TTx)]dx 

0 
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Let  1    -  ttx  =  TTy,  and  we  find 

1/tt 

f         2 
I  =      /     y      exp(-l/y)sin(l/y-7r)dy 

fl/TT    .0 

y       exp(-l/y)sin(l/y)dy 
Now  let  z  =  1/y,  and  we  have 

I  =  -  ;      e"z  sin  z  dz  =  l^e17  =  .0216070 
and  the  adequacy  of  the  evaluation  is  evident. 

We  now  consider  effective  means  of  evaluating  the  numbers  C     which 
satisfy  A-7  and  will   be  bunched  more  and  more  closely  together  as  n^°°. 
It  may  be  anticipated  that  difficulties  may  be  encountered  when  using 
customary  methods  of  root  evaluation.     In  the  applications  with  which 

this  monograph  is  concerned,  g(x)  -*  <*>  as  x    ->  b~  because  of  the  factor 

-1/2 
(1   -  x/b)         .     If  we  assume  that  in  the  interval    (a,b)  we  have 

[g(x)]2  =      Effm(x) 

where  m(x)  is   "approximately  constant,"  we  are  led  to  the  iterative 
formula 

cJ+1       bn*  (cJ-a)+a(b-cJ   )  [g(C^)]2      ^ 
n  nu  (CJn  -  a)  +  (b  -  CJn)   [g(CJn)]2  n 

where  CJ   is  the  jth  iterate  for  C   .     One  may  conveniently  take 

1  oo 

c1  =  c    ,  =  c    , 

n         n-1         n-1 


and 


c]   =  \  (a  +  b) 


However,  we  have  not  tested  this  procedure  by  actual  computations. 

For  the  case  of  two  poles  and  one  zero,  explicit  evaluations  for  Cn 
are  easily  obtained,  as  is  shown  in  Appendix  B  hereof. 
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APPENDIX  B 


The  Case  of  Two  Zeros  and  One  Pole 


Consider  the  special  case  where 

m2(s)  =  k2(s+a1)  (s+a2)/(s+b1 ) 
where  k  >  0  and  0  <  a-,  <  b,  <  a2.  We  have 


y(0)  =  yQ  =  k/  a1a2/b] 


(B-l) 


(B-2) 


and  the  inversion  is 


F(t)  =  e'yo  -  I     \       le"xtsin  y  dx  -  \    \     ^e~xtsiny    dx  (B-3) 


j       x 

ai 


IT       J  X 

a2 


In  this  case  we  are  able  to  obtain  an  explicit  representation  of  the 
values  C     introduced  in  Appendix  A.     In  both  regions  of  integration  we 


have 


y(x)  =  k/  (x-a1)(a2-x)/(b1-x) 


so  that  the  condition 


y(Cn)  =  n^ 


(cf.   equation  (A-6))  leads  to 

Cn  =  \    [(al   +  a2  +  Yn)  ~  /(a1   +  a2  +  y/  -  4(3^^^^ 


where 


Yn  =  n27r2/k 


(B-4) 


(B-5) 


(B-6) 


(B-7) 


The  lower  sign  (-)   is  to  be  used  for  the  first  integral    (from 
a]   to  b, )  and  the  upper  sign  (+)   is  to  be  used  for  the  second  integral 
(from  a2  to  °°). 
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APPENDIX  C 

Some  Additional    Inversions 

We  have  also  obtained  the  following  inversions  which  are  not 
related  to  those  considered  in  the  body  hereof  and  which  are  not  listed 
in  Reference  [  3l .     They  are  reported  here  for  record  purposes. 

f(s)  F(t) 

cosh/  as  e    cosh/  ac  2 


(s-c)   y~s  sinh/  bs  /~c"  sinh/  be 


2     y     (-1)ncos(n7r/a7b  e-n27r2t/b 
/T  ^    c  +  nV/b 


_n  a/s  (n-l)/2  

s  ne  (t/a)  I^tf/TT)         (n*l) 


Additionally,  we  note  a  rather  obvious  erratum  in  7.2.1.2,  p.  284,  of 
Reference  [3]  where  the  factor  -2s  preceding  the  sign  of  summation  should 
be  -2ir. 
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